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Abstract. - Monte Carlo investigations of magnetization versus field, M C (H), at the critical 
temperature provide direct accurate results on the critical exponent <5 _1 for one, two, three and 
four-dimensional lattices: 6^=0, <5~^=0.0666(2)~1/15, (5~^=0.1997(4)~l/5, <5~^=0.332(5)~l/3. 
This type of Monte Carlo data on 8, which is not easily found in studies of Ising lattices in 
the current literature, as far as we know, defines extremely well the numerical value of this 
exponent within very stringent limits. 



Introduction. The spin- 1/2 Ising model pQ has been a remarkably succesful model 
of a short-range interacting system to study phase transitions in magnetic, order-disorder, 
ferroelectric systems, etc. In this model, the spin variable Sj is allowed to take values ±1 
on each lattice site. The Hamiltonian that rules the interaction between the spins has the 
following form 

<ij> i 

The phase transition of the system results as a consecuence of the cooperative behavior of 
the spins determined by the first term of this hamiltonian, where J is the exchange energy. 
The second term contains the magnetic field H, which also influences the ordering of the spins. 

The critical exponents of the Ising universality class have been extensively investigated for 
a long time by various methods, High (HT) and Low Temperature (LT) expansions, Monte 
Carlo (MC) simulations, Field Theoretical (FT) methods, etc. In particular, near the critical 
point T — >T C and H — >0, the susceptibility (x) exponent (7 = <91ogx -1 /91og|T — T c |), the 
correlation length (£) exponent (y — 91og^ _1 /91og|T— T c \), the pair correlation function [r(r)] 
exponent 7/ (where T(r) / r d ~ 2+ri ) , the specific heat exponent (cv — <91ogC -1 /<91og|T — T c \), 
and the spontaneous magnetization exponent (/3 = dlogM / d\og\T — T c \) have been repeatedly 
investigated by MC methods, but no direct MC investigations of the critical isotherm exponent 
(5 = dlogH /dlogM) have been published, as far as we know. 
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Fig. 1 - Magnetization (M) vs. field (H) at T=T c —0 for a one- dimensional (d—1) Ising chain. The 
invariant magnetization confirms the expected S~ D —0 value for the critical exponent. 



Monte Carlo Method. - The determination of the critical exponent 5 from Monte Carlo 
data at the critical isotherm, M C (H) at T — T c , was made in such a way as to minimize (1) 
finite size effects, appearing only at H very close to zero, and (2) saturation effects, appearing 
at the opposite end of the sets of data, i.e. at H away from zero, when the power law M^H 1 ^ 
loses its validity. 

Relatively large lattices (of the order of 10 6 spins or more) with periodic boundary condi- 
tions and very closely spaced field intervals (AiJ~0.0005) were used to approach the critical 
point. 

In order to confirm the value of the critical temperature for each dimensionality, M s (T) was 
first determined by means of a Wolff (modified Swendsen-Wang 0) algorithm. However, 
M C (H) was studied with a standard Metropolis 0] algorithm which produces an excellent local 
and global thermalization of the whole lattice for one single temperature (in our case always 
T = T c ). In all cases, 100,000 Monte Carlo steps were used in order to insure good thermal 
equilibrium, and 20,000 states were considered in the partition function at each temperature 
/ field. 

It is easy to understand that our Metropolis algorithm, modified for studying the evolution 
of magnetization with magnetic field, is sufficient to investigate S. When we make calculations 
modifying H by very small amounts one would not expect too large changes in the fluctua- 
tions, making it unnecessary to use modified clusters algorithms, which get statistically more 
independent states for the evolution with temperature, as it is well known. In any case, the 
good quality of the results shown below, obtained with the Metropolis algorithm, speak for 
themselves. 

The M C (H) data allow the direct determination of 5 _1 from log-log plots of M vs. H at 
T = T c . 

One- dimensional Ising critical isotherm. The one-dimensional case is the easiest of 
all to be considered because the critical temperature is exactly zero (there is no real phase 
transition). In Figure 1, we show that the smallest non-zero value of the magnetic field H make 
the whole lattice turn in the same direction, perfectly ordered, keeping the magnetization at 
the maximum value (M=l) all the way. 
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Fig. 2 - Magnetization (M) vs. field (H) at T=T C =2. 269185314213 for a two-dimensional (d=2) Ising 
squared monolayer. The slope gives the 6^=0. 0666(2) value for the critical isotherm exponent. 



Two-dimensional Ising critical isotherm. - In 1944, Onsager [H] performed the calculation 
of the exact partition function of the two-dimensional Ising model in zero field. However, the 
same problem in a magnetic field remain unsolved although we can study their behavior 
numerically. In Figure 2, numerical Monte Carlo M C (H) data allow the direct determination 
of 5^=0.0666(2) from log-log plots of M vs. H at T=T C =2. 269185314213. 

Three-dimensional Ising critical isotherm. The Ising criticality for dimension d=3 
has been studied using several theoretical approaches (including MC simulations) that have 
resulted in an (almost complete) set of safe estimates of the critical exponents. Some exponents 
like 7 or f3 have received special attention and their best estimates have been summarized |B] 
as 73£>~1. 2372(5) and /3 3 £>~0. 3265(3). However, in the case of the critical exponent S^d, 
there is no direct MC measure of its value and the currently accepted estimate have been 
determined using the scaling relation 7=/3(<5-l) resulting in a value of <53r>=4. 789(2). 

As it is shown if Figure 3, our Monte Carlo M C (H) data allow the determination of 
5^5=0.1997(4) from log- log plots of M vs. H at the currently accepted value [7] of the 
critical temperature T—T C =A. 511523785. Note that this direct value of the critical exponent 
5sd—5 is appreciably different from the actual [Sj estimated value 4.789(2). 

Four- dimensional Ising critical isotherm. - It is important to note at this point, that for 
dimensions d>4 the exponents of the Ising model become the same and take the correspond- 
ing mean-field values [Sj. These exponents suddenly lock into a set of values that become 
independent of the dimensionality. 

In Figure 4, our Monte Carlo M C (H) data allow the determination of =0.332 (5) from 
log-log plots of M vs. H at the currently accepted value 9^ of the critical temperature 
T=T C =6. 68029(9). It must be noted that finite size effects begin to be observable just at 
iJ<0.001 for the particular lattice size used. 

Concluding remarks. - Table I summarizes our results for Ising systems (l<<f<4) showing 
excellent agreement between our MC results for 5 and the respective integer values 5(d) 
indicated. 
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Fig. 3 - Magnetization (M) vs. field (H) at T=T C =4.51 1523785 for a three-dimensional (d=3) Ising 
cube. The slope gives the 5^=0.1997(4) value for the critical exponent. 
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Fig. 4 - Magnetization (M) vs. field (ff) at T=T C =6. 68029 for a four-dimensional (d=4) Ising hyper- 
cube. The slope gives the 57^=0.332(5) value for the critical exponent. 



Table I - Numerical values for transition temperature (T c ) and critical isotherm exponent (8~^[ C ) for 
Ising systems with l<d<4- * Note that the value for d=3 is appreciably different from the currently 
accepted value 8sd=4-789(2) [Sj. 



d 1 2 3 4 

~f c 2.269185314213 4.511523785 6.68029 

8m C 0.0666(2) 0.1997(4) 0.332(5) 

5(d) oo 15 5* 3 
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